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ABSTRACT 

As is known, resonant scattering can distort the surface-brightness profiles of clusters 
of galaxies in X-ray lines. We demonstrate that the scattered line emission should 
be polarized and possibly detectable with near-future X-ray polarimeters. Spectrally- 
resolved mapping of a galaxy cluster in polarized X-rays could provide valuable in¬ 
dependent information on the physical conditions, in particular element abundances 
and the characteristic velocity of small-scale turbulent motions, in the intracluster gas. 
The expected degree of polarization is of the order of 10% for richest regular clusters 
(e.g. Coma) and clusters whose X-ray emission is dominated by a central cooling flow 
(e.g. Perseus or M87/Virgo). 
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1 INTRODUCTION 

Although the intergalactic gas in clusters of galaxies is op¬ 
tically thin to Thomson scattering (tt < 0.01) for the con¬ 
tinuum photons, it can be moderately thick (r ~ 1) in the 
resonance X-ray lines of highly ionized atoms of heavy ele¬ 
ments. This causes the radiative transfer in such lines to be 
important and gives rise to three major observable effects as 
previously discussed in the literature. 

First, the surface-brightness profile of a cluster in a res¬ 
onance line, calculated in the optically-thin limit, should 
be distorted due to diffusion of photons from the dense 
core into the outer regions of the cluster (Gilfanov, Sun¬ 
yaev & Churazov 1987a; Shigeyama 1998). This must be 
taken into account when attempting to determine element 
abundances from X-ray spectroscopic observations of galaxy 
clusters or haloes of elliptical galaxies. This effect has been 
already considered in application to ASCA, Beppo-SAX and 


of a bright X-ray emitting background AGN. The strongest 
absorption lines can have equivalent widths of several eV 
(Shapiro & Bahcall 1980; Basko, Komberg & Moskalenko 
1981; David 2000). By comparing the equivalent width of 
such a line and the flux from the same line seen in emission, 
the distance to the cluster can be directly determined, inde¬ 
pendent of the standard extragalactic distance scale (Krolik 
& Raymond 1988; Sarazin 1989). 

However, there should be yet another effect related to 
the problem under consideration. Namely, scattering in cer¬ 
tain (as discussed below) resonance X-ray lines causes the 
scattered emi ssion to be polarize d. The recent advent of new 
technologies (Costa et al. 2001) that promise a drastically 
increased sensitivity of near-future X-ray polarimeters has 
motivated us to consider this effect in detail. 


2 THE MODEL 


Molendi et al. 1998; |B6hringer et al. 2001; Mathews et al. 

We first consider a model cluster of galaxies which contains 

2001). 

gas with an isothermal beta-law radial density profile (Cav- 


Resonant scattering also determines the spectral profile 
(which becomes saddle-sh aped at small projecte d radii if r )> 
1) of the emergent line (Gilfanov et al. 1987a). With next 
generation X-ray spectrometers, combining high sensitivity 
with high energy resolution {E/AE ~ 10 3 TO 4 ), it should 
be possible to measure the profiles of individual lines and 
thus obtain a good deal of new information on the physical 
conditions in the intergalactic gas. 

Finally, resonant X-ray lines from the intracluster gas 
should be observable as absorption lines in the spectrum 


N e = No 1 + 


2 \ - 3 / 3/2 


(1) 


where No is the central electron number density and r c is 
the core-radius. 

At plasma temperatures typical of clusters of galaxies 
(T ~ 10' -10 8 K), all interesting X-ray lines have nearly 
Doppler absorption profiles whose width is determined by 
the velocities of thermal and bulk (turbulent) motions, since 
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the line natural width is relatively small. Indeed, the char¬ 
acteristic Doppler width is given by 

1/2 


A£V> = En 


= E 0 


2 kT 
ArripC 2 




2 kT 
Am p c 2 


(1 + 1AAM 2 ) 


1/2 


( 2 ) 


where E o is the rest energy of a given line, A is the atomic 
mass of the corresponding element, m p is the proton mass, 
Vturb is the most probable turbulent velocity (we consider 
here only small-scale motions for which the characteristic 
size of a turbulent cell is much less than the photon mean 
free path and make a simplifying assumption that the tur¬ 
bulent velocity distribution is Maxwellian), and M is the 
corresponding Mach number. Compare, for example, the ra¬ 
diative width of the Lycr line of H-like iron ( Eq = 6.97 keV), 
FWHMn = T/(2n) ~ 0.2 eV, with the corresponding 
Doppler width, 

-l 1/2 


FWHMo = 2iAn2 A E 0 = 6.0 


kT 


lOkeV 


(1 + 78M 


Therefore, the total cross section for resonant scattering 
in a given line is 


a(E) = oo exp 


/ E - Eg y 
\ ABd ) 


Here the cross section at the line center is 
v4r/ir e c/ 


do 


ABd 


(4) 


(5) 


where r e is the classical electron radius and / is the oscillator 
strength of the give n atom ic transition. 

Gilfanov et al. ( 1987a ) give a convenient expression for 
the optical depth (from the cluster center to the observer) 
at the midpoint af a Doppler-broadened resonance line: 


to 


2.7 
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where T represents the gamma-function, 5 / 8 fb ,solar is the 
abundance of a given element relative to the solar abundance 
of iron, N Z} o is the central number density of a given ion and 
iz(T) is its relative abundance at temperature T. 

The undisturbed (to -C 1) surface-brightness profile of 
line emission from the beta-cluster ([]]) is described by the 
well-known formula 


where eo = NoN Z: qA(T) [erg cm -3 s -1 ] is the volume emis- 
sivity of the plasma at the cluster center in a given line, and 
p is measured in units of r c throughout the text. 


2 \ — 3 / 3 + 1/2 


(7) 


eV.(3) 


relative weight wi) and dipole (Rayleigh) scattering (with 
W 2 = 1 — wi). The first of these components introduces no 
polarization, while the second does change the polarization 
state of the radiation field. The weights w i and w 2 are func¬ 
tions of the total angular momentum j of the ground level 
and the A j (= ±1 or 0) involved in the transition. 

The treatment of the radiative transfer problem of de¬ 
termining the surface brightness and polarization of line 
emission from the beta-cluster ([j ) is relatively simple in the 
optically-thin limit, when only the contribution of the first 
scattering is important. Although this approximation is gen¬ 
erally invalid for the cluster core region, it is always appro¬ 
priate for projected radii p > 1. We shall present below 
some analytic estimates obtained in this approximation, ad¬ 
ditionally assuming constant element abundances through¬ 
out the cluster as well as complete energy redistribution in 
scattering. The latter means that the energy distribution of 
photons after a single scattering is described by the same 
Gaussian characterizing the absorption profile (Ul). 

It follows from the spherical symmetry of the model that 
the emergent radiation will be polarized (provided that the 
given line has a dipole scattering channel) in the direction 
perpendicular to a given projected radius vector from the 
cluster center. We can therefore characterize the observed 
emission by two energy-integrated Stokes parameters: Q and 
the surface brightness B. In the single-scattering approxima¬ 
tion, the latter may be represented as B « Bo + B\ n + B 0ut , 
where Bo is the undisturbed surface brightness (Q) and Bj n 
and B ou t describe the contribution of emission singly scat¬ 
tered into and out of our line of sight, respectively. In order 
to find the Q, B lu or B au t for a given projection point p we 
need to integrate the scattered emission over all positions 
(p, l ) along the line of sight to the observer (here l is the 
coordinate along the line of sight). 

We thus may write for isotropic scattering 

Bin(p) = J dlN z j dc °*° d * l'(p,l,e,<j>), Q( P ) = 0.(8) 

For dipole scattering we have 

Bin{p) = -^5= J dlN z J d cos 6 dct> 3(1 + 6 C ° s2 9) I'{p , l, 9, (A)(9) 
and 

Q(p) = [ dlN z [ d cos 6 d(j> ^ cos(20)7 / (p, Z, 9, 0). (10) 

V2 J J 

In the above formulae the factor 1 /\[2 results from the con¬ 
volution of the incident (Gaussian) line energy profile with 
the energy-dependent cross-section (Q) , and I' is the solution 
of the undisturbed problem for energy-integrated intensity 
(see below). We have not given above the corresponding for¬ 
mulae for Bout, which are similar to equations (|s|) and (|)|) 
for Bin, because we shall not use them explicitly. 


2.2 Calculations 


2.1 Scattering phase function and polarization 

According to quantum mechanics, resonant scatteri ng can 
be represented as a combinati on of two processes (Hamil- 
ton 1947; Chandrasekhar 1950): isotropic scattering (with a 


Below we present some results obtained in the single¬ 
scattering approximation using equations (|i]) f p~C)|) an d ap¬ 
pend some resu lts obtained e arlier by Sunyaev (1982) and 
Gilfanov et al. ( 1987a ; 1987b ). We then compare these ana¬ 
lytic expressions with the results of numerical computations 
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that allow for multiple resonant scatterings which become 
of importance when to > 1. 

Our Monte-Carlo code allows one to treat an individ¬ 
ual act of resonant scattering in full detail, i.e. not re¬ 
sorting to the hypothesis of complete energy redistribu¬ 
tion. This is realized in a straightforward way: for a pho¬ 
ton propagating in a given direction fl with given energy E 
and polarization direction m, first an ion is drawn from a 
Maxwellian velocity distribution which has a velocity com¬ 
ponent (vfi) = c(E — Eq)/Eo (in the limit v <S c). Then 
the direction fi' of the emergent photon is drawn in accor¬ 
dance with the relevant scattering phase matrix (isotropic 
or dipole), and the corresponding (Doppler-shifted) energy 
E' = E 0 [l + (vfi')/c] is found. 


2.2.1 Ideal beta-cluster 


Isotropic scattering 


Gilfanov et al. ( 1987a ) have derived the surface-brightness 
profile of the beta-cluster (lj) in the limit to < 1 for par¬ 
ticular (3 values and explicitly assuming an isotropic phase 
function: 


B(J3 = 2/3) = 


r c eo 


1 


to 


To 


(1 + p 2 ) 3 / 2 \[2 (1 + p 2 ) 2 2^2(1 + p 2 ) 3/2 


( 11 ) 


m= i) 


3r c e 0 

32 


1 r 0 1 V2t 0 13 + 5p 2 

(1+p 2 ) 5 / 2 “ 7^(1 +P 2 ) 7/2 + “45“ (1 + p 2 ) 3 


( 12 ) 


Note that the above expressions take into account the (neg¬ 
ative) contribution of emission scattered from our line of 
sight into other directions. 


Dipole scattering 

In this case we must substitute into equations (|)j) and ( |To[ ) 
the solution of the undisturbed problem, 


7(r„,a), 


ds 


, r 2 
1 3—2 
r£ 


- 3/3 


(13) 


where the integration is done along the direction fl = 
(sin 0 cos <j>, sin 9 sin <j>, cos 6) and r = ro + sEl. 

It proves impossible to do analytically the final inte¬ 
gral (over the line of sight) in equations ([)]) and (R). How¬ 
ever, the integration can be completed in the limit /) > 1 
if (3 > 1/2. The latter condition, which means that the to¬ 
tal luminosity of the cluster converges when its outer radius 
r ou t —> oo, ensures that for a scattering site located far (at 
r ^$> r c ) from the cluster center the flux from the core region 
subtending a small solid angle is much larger than the cumu¬ 
lative flux of radiation coming from all othe r directions. Note 
that (3 > 1/2 for most observed clusters (lones & Forman 


199E). The resulting asymptotes for p max(l,r 0 




3 f2 r(3/3/2)T(3/3 - 3/2) 

64 V t V(3(3/2 - l/2)r(3/3/2 + 2)T(3/3) 


x [T(3/3/2 + 1/2) + F(3/3/2 + 3/2)] r c e 0 r 0 p“ 3/3 - 1 , (14) 


0 = 


3 / 2 T(3/3/2)T(3/3/2 + 3/2)T(3/3 — 3/2) 


64 V tt r(3/3/2 - l/2)F(3/3/2 + 2)T(3 /3) 


xr c eorop 


- 3 / 3-1 


1 T(3/3 — 1/2) - 6 / 3+1 

Bo ~WH rm C0P ■ 

The degree of polarization can then be found as 

Q _ Q 


P = 


Bo + Bin + Bout Bo + Bj, 


(15) 

(16) 


(17) 


Note that the above formulae are applicable (for suf¬ 
ficiently large p) for arbitrary values of the central opti¬ 
cal depth ro, not only in the limit tq < 1. Indeed, outer 
cluster regions are optically thin to resonant scattering 
(r(p) ~ Top~ 3f3+1 ) and the amount of emission scattered 
there is simply proportional (for j3 > 1/2) to the total clus¬ 
ter luminosity dominated by emission from the central dense 
core, regardless of whether this core is optically thin or thick. 
For this reason it is important to include Bi n in the de¬ 
nominator of the expression ([t7j), as its contribution to the 
observable surface brightness can be large or even domi¬ 
nant. On the other hand, the corresponding contribution 
from emission scattered out of our line of sight is relatively 
small: B out ~ r(p)B 0 (p) < B in when (3 > 1/2. 

In particular, we find from equations ( [ill ) (0) that 

P( 3 = 2/3) = (9^2/64)t 0 0.20r o 

1 + (15\/2/64)to 1 + 0.33ro 

for p 3> max(l, ro) (18) 


and 

P(/3 = 1) 


(V2/15)r 0 p _ 0.094rop 

1 + (v/2/10)r 0 p ~ 1 + 0.141r o p 

for p max(l, t^ 2 ). 


(19) 


The polarization quickly vanishes when one approaches 
the cluster center (p —> 0). This statement is readily proved 
for ro 1, when we have 


P(p< 1,(3 = 2/3) 


9(3-4 In 2) 

32^2 


2 


Top 


0.045r 0 p 2 


( 20 ) 


and 


/2 

P(p < 1,0 ^ 1) = -—Top 2 w 0.094rop 2 . (21) 

15 

We now proceed to discussing the results of our nu¬ 
merical computations. It turns out that the cluster surface- 
brightness profile (not the polarization!) corresponding to 
a given set of parameter values is almost insensitive to the 
scattering phase function used, isotropic or dipole. For this 
reason we shall not present below any radial profiles cor¬ 
responding to isotropic scattering and concentrate on the 
dipole scattering case. 

Figures ([!]) (]3|) show a set of Monte-Carlo radial pro¬ 
files of the Stokes parameter Q and surface brightness B as 
well as the degree of polarization P for two (3 values: 2/3 
and 1. We see that P grows monotonically with increasing 
p, being negligible near the cluster center and approaching 
at o 1 the asymptote given by equation ( fl8| ) or equation 
(|19|) if /3 = 2/3 or /3 = 1, respectively. In the former case, this 
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Figure 1. Degree of polarization in a resonance line as a func¬ 
tion of projected distance from the center of a beta-cluster with 
(3 = 2/3 for different optical depths at the line center (the solid 
curves). The dashed curves are the corresponding p —> oo asymp¬ 
totes calculated from equation (|l8f)- 



Projected radius 


Figure 2. Total surface brightness ( B ) and its polarized compo¬ 
nent (Q) in a resonance line as functions of projected radius for 
a beta-cluster with (3 = 2/3 for two values of the optical depth at 
the line center. The dotted line shows the undisturbed (to = 0) 
surface brightness profile. 



Figure 3. Same as Fig. |l|, but for (3 = 1. The dashed curves are 
p —> oo asymptotes calculated from equation (|19[). 



Projected radius 


Figure 4. Same as Fig. ^ but for (3=1. 

asymptote is a constant, while when (3 = 1, the polarization 
grows proportionally to p until the projected distance be¬ 
comes sufficiently large that the radiation scattered into our 
line of sight begins to dominate the surface brightness (see 
the denominator of the expression (|l^)). The polarization at 
a given p also grows with increasing optical depth to until, 
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Projected radius 


Figure 5. Degree of polarization in a resonance line as a function 
of projected distance from the center of a beta-cluster for different 
values of the /3-parameter and ro = 1. 

again, at a certain value of ro ~ 10 the scattered emission 
becomes prevailing and P does not increase anymore. 

The effect of multiple resonant scatterings is apparent 
in the central region of the cluster (p < a few), where the 
degree of polarization decreases as p —> 0 even faster than 
given by the asymptotic formulae (^(J) and (j2l]) when to > 1 . 

Let us now explore how the polarization effect depends 
on the /3-parameter when its value is allowed to vary in a 
broader range. Fig. compares a set of radial profiles of P 
obtained for different /3 values and a given ro = 1. As ex¬ 
pected, the degree of polarization diminishes with decreasing 
/3 as the cluster emission becomes less centrally peaked and 
the radiation field permeating the intracluster gas becomes 
more isotropic. Far from the cluster core, the degree of polar¬ 
ization is a weak declining function of p when 1/2 < fl < 2/3. 
Indeed, in this case P(p 33> 1) ~ Q/Bo ~ rop 3#3 ~ 2 in accor¬ 
dance with equations 0 5) (|l7|). For example, P{p l,/3 = 
0 .6) ~ r 0 p-°' 2 . 

The degree of polarization remains a declining function 
of p 1 also when /3 < 1/2. In this case the core region 
no longer dominates the radiation field within the cluster 
and polarization is mainly caused by the large-scale angular 
anisotropy of the radiation field. Although this case is not 
easily analytically treatable, we may expect that Q ~ Bi n ~ 
t(p)B 0 and P(p » 1) fts Q/B 0 ~ r(p) ~ r 0 p~ 3/3+1 . For 
example, P(p S> 1, /3 = 0.4) ~ rop~°' 2 , the same dependence 
as in the case of /3 = 0.6. It is clear from the profiles shown 
in Fig. [i| that when /3 < 2/3, the polarization degree can 
be roughly described by a constant value at p j> 5 which 
depends on /3 and ro. 

Given the fact that the degree of polarization is an in¬ 
creasing function of p at least within a few core radii from 



0.5 1 5 10 


Projected radius 


Figure 6. Signal-to-noise ratio (defined as Q/pl3) in detecting 
polarization in a resonance line from a beta-cluster as a function 
of projected radius for two valus of the central optical depth and 
different /3 values: 2/3 (the solid line), 1 (the dashed line) and 0.5 
(the dotted line). 

the cluster center while the surface brightness falls off with 
moving away from the cluster core, we may anticipate that 
there should exist a range of projected radii favorable from 
the point of view of detectibility of the polarized signal. In 
order to check this suggestion we have plotted in Fig. [] the 
signal-to-noise ratio, which is proportional to Q/\/B, as a 
function of p. In doing so we, of course, assumed that the 
contribution of the instrumental and cosmic background to 
the detected X-ray flux is negligible. One can see that the 
region 0.5 <( p <( 10 centered (on a logarithmic scale) at 
po ~ 2 is most suitable for measurements with future X-ray 
polarimeters, with the expected degree of polarization being 
~ 10/3 2 to per cent at po = 2 for 0.5 < /3 < 1 and ro < 2. 

Although in obtaining the numerical results described 
above we treated the process of resonant scattering accu¬ 
rately, assuming complete energy redistribution in scatter- 
img leads to virtually the same results. This is demonstrated 
in Fig. [i], where we have plotted the degree of polarization 
as a function of photon energy for different projected radii. 
Note that in the considered case of a moderate optical depth 
(ro = 4), P is somewhat higher in the wings of the resonance 
line than at its center for the spectrum taken from the cluster 
core region, while the polarization is approximately constant 
across the line spectra measured at p 1. 

2.2.2 Beta-cluster with a central cooling flow or AGN 

We might expect that the net degree of polarization would 
be larger for clusters with a luminous cooling flow. Observed 
X-ray surface brightness profiles of cooling flows usually im¬ 
ply an approximately 1/r radial density law for the X-ray 
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(E-E 0 )/AE„ 


Figure 7. Degree of polarization as a function of photon energy 
across the spectrum of a resonance line with To = 4 for a beta- 
cluster with 3 = 2/3 at different projected distances from the 
cluster center. The solid lines result from accurate Monte-Carlo 
calculations. The dashed lines were obtained assuming complete 
energy redistribution in scattering. The curves corresponding to 
the innermost region p < 0.5 were multiplied by 4 to make the 
trough centered on E = Eq visible. 


emitting gas which probably flattens out at r < a few kpc. 
Therefore, when considering a cooling-flow beta-cluster, one 
possibility is to represent the cooling flow by an additional 
beta-model with 0 rs 1/3 and a small core-radius r c . 

At the same time, it is illuminating to consider this situ¬ 
ation as though we have a beta-cluster with total luminosity 
Lp in a given line and a point source in the cluster center 
with a luminosity L in the same line. Clearly, the results 
of such modelling will be applicable only at projected radii 
larger than the effective size of the cooling flow (typically at 
p > a few tens of kpc). 

At the same time the results reported below are applica¬ 
ble to another important situation in which there is a strong 
active galaxy (AGN) in the center of a beta-cluster. In this 
case L will represent the luminosity of resonantly absorbed 
radiation from the AGN continuum (power-law) spectrum 
in a given X-ray line. 

To find the contribution of the central point source 
to the surface brightness and polarization in the single¬ 
scattering approximation we must substitute 

J(r 0 , f2) = -t— 2<5(fi — ro/ro), (22) 

47T 7q 

into equations (jij) (^o|), which leads to the results given be¬ 
low. 


Isotropic scattering 
If f3 = 2/3, 


= ( 1 

m 167T 2 r 2 ip 

If 0 = 1, 

D _ V2 -I/To ( 7T 

in ~ 32 tt 2 r 2 \p 



Q = 0. 


(23) 


2 2arctanp\ 

TT^ p ) 


<2 = 0. (24) 


Dipole scattering 
If 0 = 2/3, 


_ 3^2 Ltq ( 3 + 2 p 2 _ 4 + 2 p 2 ^ 
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B h 


3\/2 Ltq 
128tt 2 r 2 


( 5 + 3 p 2 
l 1 +P 2 


P P 


Q 


3V2 Lt 0 
1287t 2 r| 


( 1 + 3p 2 

V i + P 2 


+ 


1 - 3p 2 1 

-arctan — 

P P 


(25) 


(26) 


One can see that the effect of a particular phase func¬ 
tion (isotropic or dipole) on the surface-brightness profile of 
scattered emission is quite small even in the case under con¬ 
sideration, where we are dealing with a compact (not diffuse 
as in the case of a self-illuminating beta-cluster) illuminat¬ 
ing source. In particular if /3 = 2/3, /Iscat = 9/8 for 
p -► 0 and 15/16 for p -> 00 . If 0 = 1, J^/4 St = 9/8 for 
p —> 0 and 9/10 for p —> oo. 

If we ignore for the moment the resonant scattering of 
the beta-cluster emission, then the total observed surface 
brightness will be B = Bo + B ln and the corresponding po¬ 
larization degree P = Q/(Bo + -Bin), with Bo being given by 
equation (jr). In the vicinity of the central source (at p< 1) 
the surface brightness will be dominated by the scattered 
emission, so that P = Q/Bi n . In the case of dipole scatter¬ 
ing, P —> 1/3 when p —> 0 for arbitrary values of 0. The 
polarization at large projected radii (p 1), will depend 
on the ratio B ln /Bo ~ (L/Lp)rop 313 2 ■ This ratio is grow¬ 
ing with distance from the cluster center for 0 > 2/3, and 
declining in the opposite case. 

If the condition B lu 3> -Bo is indeed met, then 


B(p» 1) = 


1 

1 + T(30/2 + l/2)/T(30/2 + 3/2)' 


(27) 


This relation follows directly from equations and 

(0 , which were originally derived for the ideal beta-cluster 
assuming that its entire luminosity was emitted from its 
center. In particular, P(p 1,0 = 2/3) = 3/5, P(p » 
1,0 = 1) = 2/3. 

Figs. H and I show the results of our numerical simu¬ 
lations for the model of a beta-cluster with a point central 
source. One can see that deep within the cluster core the 
degree of polarization approaches the value 1/3 as p —> 0, 
in agreement with the analytic result, because the radiation 
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field is completely dominated by emission from the point 
source. When we move away from the core region to p ^ a 
few, provided that (3 > 1/2, we effectively come to a sit¬ 
uation in which there is a point central source with total 
luminosity Ltot = L + L f 3 whose emission is scattered by the 
beta-cluster. The resulting P(p) profile is then a function of 
L/Lq. If (3 > 2/3, P approaches as p —> 00 the asymptotic 
value characterizing the point-source case: 3/5 if (3 = 2/3 
and 2/3 when /3 = 1. If (3 < 1/2, the role of the central 
point source becomes unimportant at p 1, si milar ly to 
the role of the beta-cluster core, as discussed in §2 .2.1 . 

We see from Figs. and || that already at L = 0.5Lp the 
observed line emission is strongly polarized across the whole 
cluster, with P being of the order of 15% when tq = 1. 


3 IMPLICATIONS FOR OBSERVATIONS 

Let us now discuss the prospects of observing the resonant 
scattering polarization effect in actual clusters of galaxies. 
We shall begin by discussing potentially interesting X-ray 
lines and then continue by performing numerical simulations 
for a few clusters. 


Figure 8. Polarization degree in a resonance line with ro = 1 as 
a function of projected radius for a j3 = 0.5 beta-cluster with a 
point source in the center for different ratios of the luminosity of 
the central source (L) to the total luminosity of the cluster Lp in 
the line. 


3.1 Most promising lines 

For the discussion below as well as for the calculations in 


i.2 it is important that given the very low density of the 


intracluster gas, virtually all ions are in their ground atomic 
state. In particular, only the lowest fine-structure sublevels 
are populated. 



Figure 9. Same as Fig. ^j, but for 3 = 2/3. Also shown is the 
analytic result (E5|) valid in the small optical depth limit (ro <C 1). 


3.1.1 He-like K a line 

The resonant transition Is 2 — ls2p( 1 Pi) at 6.7 keV of He-like 
iron corresponds to one of the most prominent lines in the 
spectra of rich clusters of galaxies. He-like iron is abundant 
in a wide range of temperatures from 1.5 to 10 keV. Taking 
into account that the transition has a large absorption os¬ 
cillator strength of ~ 0.7, the expected optical depth in the 
line can be large, and therefore this line is one of the prime 
candidates for observing the resonant scattering effect. Since 
the transition corresponds to the change of the total angular 
momentum from 0 to 1, the phase function of the scattering 
is a pure dipole. This means that this line is also expected 
to have highest degree of polarization. 

In the vicinity of this line there are several other lines 
including the intercombination Is 2 —ls2p( 3 P) and forbidden 
Is 2 — ls2p( 3 S) lines of He-like iron and many dielectronic 
satellite lines. For the forbidden line the optical depth is 
essentially zero, while for the most important intercombina¬ 
tion line (Is 2 — ls2p( 3 A)) the absorption oscillator strength 
is at least a factor of 10 lower than for the resonant line. 
While many of the satellite lines have a vanishing optical 
depth, some satellite lines do have an appreciable oscillator 
strength. The latter correspond to the singly exited state of 
particular ions of iron (e.g. the transitions ls 2 2s — Is2s2p 
for Li-like iron). However, Li-like iron is present in sufficient 
amounts only in a relatively cool (below 4 keV) plasma, 
which implies that in hot clusters the optical depth in these 
lines will be lower than for the resonant He line. Note also 
that for the satellite lines autoionization of electrons (from 
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Figure 10. Decrease of the observed degree of polarization for 
an instrument collecting photons in a given energy range centered 
on the He-like line at 6.7 keV. The solid and dotted lines corre¬ 
spond to a plasma temperature of 8 and 2 keV, respectively. The 
abundance of iron was assumed to be 0.5 solar. Only the flux in 
the line is assumed to be polarized. The edges in the curves cor¬ 
respond to such widths of the energy window when an additional 
prominent line falls inside the window. 

the exited state) plays a role, thus suppressing the emis¬ 
sion in the satellite line in the case of a sufficient optical 
depth. A similar effect called “resonant Auger destruction” 
has been considered in application to ionized accretion disks 
(e.g. Ross, Fabian & Brandt 1996). For clusters this effect 
is less important because of the relatively low optical depth 
and lower ratio of autoionization and radiative decay rates 
(at least for relevant transitions in the Li-like iron). 

All these lines (except for dielectronic satellites with 
a large principal number of the spectator electron in Li- 
like iron) are separated from the main resonant line by at 
least 30 eV. Therefore, even turbulent broadening does not 
severely mix the resonant line with other lines. Finite energy 
resolution of the polarimeter may, however, significantly af¬ 
fect the prospects of detection of the polarized emission. 
Assuming that only the resonant line of He-like iron is polar¬ 
ized, the measurable polarization can be estimated by multi¬ 
plying the degree of polarization as would result in the ideal 
case (as calculated in §'/) by the factor / = F\- me /Ft„t where 
Fline is the flux in the resonant line and Ltot is the total flux 
measured by the detector in an energy band A E centered 
on the resonant line energy. In Fig. [Rj we plot this factor as 
a function of A E for two values of the plasma temperature 
of 8 and 2 keV. One can see that the energy resolution of a 
CCD-type detector is needed to avoid a strong decrease in 
polarization due to contamination by the unpolarized flux 
in the continuum and other lines. 


3.1.2 Fte-like Kp iron line 

This line (Is 2 — ls3p( 1 Pi)) at 7.88 keV is very similar to 
the K a line, but the oscillator strength is ~ 0.16, i.e. a fac¬ 
tor of 4.5 smaller. This directly translates into the optical 
depth being smaller by the same factor. It was suggested 
(e.g. Akimoto et al. 1997; Molendi 1998) that the compar¬ 
ison of the radial profiles in the K a and Kp lines may be 
used as an indicator of the resonant scattering in clusters. 
Note, however, that the different energies of these two lines 
imply some dependence of the ratio of their emissivities on 
temperature. This effect may complicate an analysis of real 
data. 

3.1.3 H-like K a iron line 

This line at 6.96 keV has two components, Is — 2p( 1 L > 1 / 2 ) 
and Is — 2p( 1 P 3 / 2 ), separated by ~ 20 eV. The oscillator 
strengths of the lines are 0.14 and 0.28 respectively. For the 
first line the total angular momentum is equal to 1/2 both 
for the ground and excited states, and therefore the scat¬ 
tering phase function is isotropic and scattered emission is 
unpolarized. For the second line the total angular momen¬ 
tum of the exited state is 3/2 and the phase function is 
a sum of isotropic and dipole phase functions with equal 
weights. For the conditions in clusters (i.e. only marginally 
optically thick cases) all complications associated with the 
quantum interference of these two lines (e.g. Stenflo 1980) 
can be safely neglected. We note here that the ratio of the 
intensities of these two lines is essentially independent of 
temperature. Therefore, a comparison of the radial profiles 
and polarization properties of these lines would be the clean¬ 
est test for the presence of resonant scattering, should next 
generation calorimeters and polarimeters allow one to re¬ 
solve them. However hydrogen-like iron is abundant only in 
high temperature clusters, and the optical depth is system¬ 
atically lower than for the He-like resonant line. 

3.1.4 K-shell lines of lighter elements 

For low temperature clusters and especially for cooling flow 
regions with the gas temperatures in the range of 1-3 keV, 
the K-shell lines of lighter elements (e.g. Mg, Si or S) may 
be important. These elements have abundances comparable 
to the iron abundance and the main limiting factor is the 
ionization fraction. Compared to iron K-lines, the cross sec¬ 
tions (per ion) for these lines are higher - see equations (^), 
(|[2!) because of the lower energy of the transition, although 
the smaller atomic masses imply that Doppler broadening 
is more important for these elements than for iron (and as a 
consequence turbulent broadening has less impact on these 
lines than on iron K-lines). 

3.1.5 L-shell lines of iron 

The L-shell transitions (e.g. of the 2s — 3 p or 2 p — 3 d type) 
with energies of the order of 1 keV are especially strong 
at lower temperatures of 1 -3 keV typical for cooling flow 
regions. They are associated with Li-, Be-, B-, C-, N-, 0-, F- 
and Ne-like ions of iron. At any given temperature several 
ionization stages are present and several strong lines are 
present. Polarization properties of these transitions depend 
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on the total angular momenta of the ground and excited 
states. 

Several lines appear particularly promising and inter¬ 
esting for polarimetric studies. At temperatures of the order 
of 1 keV this is the Be-like iron line at 1.129 keV which 
has a large oscillator strength (/ = 0.7) and a pure dipole 
scattering phase function. At still lower temperatures ~ 0.5 
keV the most interesting is the Ne-like iron (Fe XVII) line 
at 0.825 keV which has an exceptionally large oscillator 
strength (/ = 2.95) and a pure dipole scattering phase func¬ 
tion. We note however that it is not clear if the gas in this 
temperature range is abundant in cluster cooling flows. The 
standard models do predict such low-tempera ture plasma, 
while recent Chandra and XMM observations ( Bohringer et 
al. 200^) do not provide clear evidence for it. 


Table 1. X-ray lines with significant optical depth (to > 0.05) to 
resonant scattering. The Coma cluster. 


Ion 

Energy 

(keV) 

Optical depth 

Weight of dipole 
scattering 

Fe XXIV 

1.163 

0.04 

0 

Fe XXIV 

1.168 

0.07 

0.5 

Fe XXV 

6.700 

0.37 

1 

Fe XXV 

7.881 

0.06 

1 

Fe XXVI 

6.952 

0.05 

0 

Fe XXVI 

6.973 

0.09 

0.5 


3.2 Individual clusters 


We have performed Monte-Carlo calculations to simulate 
the resonant scattering polarization effect in three proto¬ 
typical clusters of galaxies. The line energies and oscillator 
strengths were taken from th e list of strong resonant lines 
of Verner, Verner & Ferland ( 1996 ). The solar a bund ances 
of elements were taken from Anders fe Grevesse (198S). We 


used the MEKAL code (Mewe et al. 1981; Mewe et al. 1986 


Kaastra 1992 ; Licdahl et al. 1995 ) as implem ented in the 
software package XSPEC vlO ( Arnaud 199t| ) to calculate 
plasma emissivity in a given line as a function of radius. The 
ionization fractions were cal culate d using collisional ioniza¬ 
tion rate fits from Voronov ( 1997 ), radiative recombination 
rates from Verner & Ferland (1996) and dielect ronic recom¬ 
bination rates fro m Ald rovandi & Pequignot (1973), Shu ll 
& van Steenberg ( 1982| ) and Arnaud & Rothenflug (1985). 
We note here that this usage of several independent sources 
of atomic data may lead to certain inconsistencies in the 
resulting values of the polarization degree, but we believe 
that the net effect (of these inconsistencies) is very minor 
and can be neglected. 


Coma Cluster 



Projected radius (kpc) 


3.2.1 Coma cluster 


The Coma cluster is the prototype of a rich regular cluster. 
The temperature of the gas is kT = 8.1 keV and its density 
distribution is described by the beta-law (Hi with (3 = 0.67, 
r c = 430 kpc and No = 2.3 x 10 -3 cm -3 (Jones & Forman 


1991). According to recent XMM-Newton observat ions the 


abundance of iron is 0.25 relative to the solar value (Arnaud 


et al. 2001). At the given gas temperature He-like and H-like 


ions of iron are present in nearly equal amounts with a neg¬ 
ligible admixture of less ionized atoms. Given the oscillator 
strengths of the corresponding lines, only the 6.7 keV line 
of He-like iron proves to have a substantial optical depth to 
resonant scattering, as indicated in Table |l|. 

The reported values of to were calculated assuming 
A4 = 0. The optical depth will be smaller is we adopt a 
more realistic M > 0 value, because turbulent broaden¬ 
ing is especially important for iron, as follows from equa¬ 
tion (|5|). For example, to = 0.18 for the 6.7 keV line when 
M = 0.2. Fig. [TT] shows the expected radial profile of the 
degree of polarization in the 6.7 keV line as would be mea¬ 
sured by a polarimeter with perfect spectral resolution. We 


Figure 11. Polarization degree in the most prominent resonant 
line as a function of projected distance from the center of the 
Coma cluster. The solid line corresponds to the case of negligible 
turbulence in the gas, while the ratio of the characteristic turbu¬ 
lent velocity to the sound velocity is J\4 = 0.2 and M = 0.4 for 
the dotted and dashed line, respectively. 


see that P < 5% is expected. The polarization would be a 
little smaller when measured by a detector having a finite 
spectral resolution, as is clear from Fig. [1^. 

3.2.2 Perseus cluster (A426) 

The X-ray emission of the Perseus cluster is peaked on the 
central dominant galaxy NGC 1275, suggesting the presence 
of a cooling flow. 

We adopted simple analytic expressions for the radial 
dependence of density and temperature in this cluster. The 
overall density distribution, including the central cooling 
flow region, was described by a single beta-model (Q) with 
the parameters No = 0.03 cm -3 , /? = 0.47 and r c = 47 kpc. 
The temperature distribution was set to 
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Table 2. X-ray lines with large optical depth (to > 0.5) to reso¬ 
nant scattering. The Perseus cluster. 


Ion 

Energy 

(keV) 

Optical depth 

Weight of dipole 
scattering 

Fe XXIV 

1.163 

0.8 

0 

Fe XXIV 

1.168 

1.6 

0.5 

Fe XXV 

6.700 

3.3 

1 

Fe XXV 

7.881 

0.5 

1 


kT = 7 


1 + {r/rpf 
2.3 + (r/r o) 3 


keV, 


(28) 


with r o = 100 kpc. The gas is coldest at the cluster cen¬ 
ter, where kT = 3 keV, and hottest in the outermost re¬ 
gions, reaching kT = 7 keV. The above distributions ap¬ 
proximately correspond to the deprojected density and tem¬ 
perature profiles derived from XMM-Newton data (to be 
published elsewhere). Given the substantial azimuthal asym¬ 
metry of the Perseus cluster, these azimuthly averaged dis¬ 
tributions can be used only for crude estimates. Note that 
outside the cooling flow (at r > 300 kpc) the density ra¬ 
dial profile we use approximately matches the beta-profil e 
inferred from Einstein observations (Jones & Forman 1999), 
defined by values Nq = 4.1 x 10~ 3 cm -1 , (3 = 0.58 and 
r c = 280 kpc. The ratio of the two profiles is less than 1.5 
up to a distance of 4 Mpc. 

We have adopted a constant abundance of iron across 
the cluster of 0.5 solar, as inferred from the XMM-Newton 
spectroscopic data. 

Table lists resonant lines for which to > 0.5 (at 

We s oo that tiro 6 . 7 koV line i s again the one with 


M = 0 

the largest optical depth Tn contrast to the previous wam- 

ple, Li-like iron is fairly abundant in the Perseus cluster, 
particularly in its cooling flow, which gives rise to a mod¬ 
erate optical depth in the Fe XXIV doublet at 1.17 keV. 
Fig. |l2| shows the simulated radial profiles for the degree 
of polarization in the lines listed in Table In the case of 
the 1.17 keV doublet, the polarization degree was found as 
P = (Qi+Q 2 )/{Bi + B 2 ) (with Qi ,2 and B\^ being the cor¬ 
responding quantities for each of the two lines). This calcu¬ 
lation assumes that the polarimeter used cannot resolve the 
lines. Note that Doppler broadening is not sufficient to cause 
these lines to overlap even at sonic turbulence (A4 = 1): 


AE d = 0.5 


kT 

5kW< 1 + 78M 


1/2 


eV, 


(29) 


compared with the 4.6 eV splitting of the doublet. 

We see from Fig. |l2j that the expected degree of po¬ 
larization is high for the 6.7 keV line - larger than 8% at 
p > 300 kpc. The polarization is a factor of 5 smaller for the 
7.88 keV and 1.17 keV lines. This is caused by the smaller 
optical depth and in the latter case also by the lower weight 
of dipole scattering. We also see that the degree of polar¬ 
ization is nearly constant outside the central 300 kpc for all 
three sampled lines, which is in agreement with the results 

of §| 


Perseus Cluster 



Figure 12. Same as 



but for the Perseus cluster. 


3.2.3 M87/Virgo cluster 


This is another example of a cooling flow cluster, but in 
this case the characteristic temperature of the cooling flow 
is lower (kT ~ 1.5 keV in the center) than in the proceed¬ 
ing example. For our computations we have adopted gas 
distribution paramete rs as inferred from XMM-Newton and 
ROSAT observations (Nulsen fe Boehringer 1995; Bohringer 


et al. 2001; Mathews et al. 2001). Namely, the density dis¬ 


tribution within the central 12 kpc is represented by a beta- 
model with (3 = 0.47, No = 0.35 cm -3 and r c = 2.0 kpc. At 
larger distances, a less steep density law is used: 


N e = 0.12 


1 + 


4.3 kpc 


(30) 


which is approximately a beta-law with fa = 0.39. 

The overall gas temperature distribution is approxi¬ 
mated by 


T = To 


1 + (r/ 13.2 kpc) 


+ 0.20 


13.2 kpc 


(31) 


The temperature rises from the center, where To = 1.3 keV, 
towards a broad maximum of T = 3.0 keV at r ~ 150 kpc. 

The element abundances are assumed to be constant 
at r < 10 kpc, being F( O) = 0.4, F(Si) = 1.0, F( S) = 1.0, 
T(Ar) = 1.0, F( Ca) = 1.0 and T(Fe) = 0.7 (in units of solar 
abundances), then gradually falling with radius to reach the 
values F( O) = 0.4 (unchanged), T(Si) = 0.6, F( S) = 0.6, 
T(Ar) = 0.6, F( Ca) = 0.6 and T(Fe) = 0.35 at r = 40 kpc, 
and remain constant from there on. 

Table |i| gives a list of potentially interesting lines for 
observing the resonant scattering polarization effect. There 
are in total 17 lines (if counting fine-structure multiplets 
as a single line) with optical depth to resonant scattering 
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Table 3. X-ray lines with large optical depth (to > 0.5) to reso¬ 
nant scattering. The M87/Virgo cooling flow. 


Ion 

Energy 

(keV) 

Optical depth 

Weight of dipole 
scattering 

O VIII 

0.654 

0.4 

0 

O VIII 

0.654 

0.7 

0.5 

Ne X 

1.021 

0.5 

0 

Ne X 

1.022 

1.0 

0.5 

Si XIII 

1.865 

2.1 

1 

Si XIV 

2.004 

0.8 

0 

Si XIV 

2.006 

1.7 

0.5 

S XV 

2.461 

2.5 

1 

S XVI 

2.620 

0.3 

0 

S XVI 

2.623 

0.6 

0.5 

Ar XVII 

3.140 

0.8 

1 

Fe XX 

0.967 

0.7 

0.28 

Fe XX 

0.967 

0.5 

0.32 

Fe XXI 

1.009 

2.3 

1 

Fe XXII 

1.053 

3.3 

0.5 

Fe XXII 

1.064 

0.5 

0 

Fe XXII 

1.084 

1.5 

0.5 

Fe XXIII 

1.129 

6.2 

1 

Fe XXIII 

1.493 

1.2 

1 

Fe XXIV 

1.163 

1.6 

0 

Fe XXIV 

1.168 

3.2 

0.5 

Fe XXIV 

1.553 

0.6 

0.5 

Fe XXV 

6.700 

1.8 

1 


exceeding 0.5. Most of these lines, including the lines of el¬ 
ements lighter than Si, the lines of He-like Si and S, and 
those of Fe ionic species less ionized than Fe XXIV origi¬ 
nate preferentially within a few tens of kpc from the M87 
center where the temperature is sufficiently low. The most 
interesting line is the one at 1.129 keV, which corresponds 
to the 2s — 3p transition in Be-like iron. This line has the 
largest optical depth and a dipole scattering phase function. 
Also notable are the Ka lines of helium-like Si and S at 
1.865 and 2.461 keV, respectively. These lines are well sepa¬ 
rated on the spectrum from the Fe L-shell complex of lines 
and the detection of polarization in them does not require a 
high spectral resolution of the polarimeter. 

Fig. i shows the simulated radial profiles of the degree 
of polarization for several prominent lines. As expected, the 
highest polarization is achieved in the Fe XXIII 1.129 keV 
line, reaching a maximum of 14% (at A4 = 0) very close 
to the cooling flow center (at p ~ 15 kpc), where the sur¬ 
face brightness is still large. The polarization degree falls off 
from this maximum with moving away from the center, and 
a similar behaviour pertains to the other lines, excluding the 
lines of Li-like (e.g. the one at 1.17 keV) and He-like iron 
(6.7 keV), which are characterized by increasing P with pro¬ 
jected radius. The case of the 1.129 keV line and other sim¬ 
ilar lines closely corresponds to our model of a poin t source 
in the center of a beta-cluster considered in §2.2.2. Indeed, 
the fraction of the overall line luminosity emitted within the 
central 10 kpc is as large as 33% for the Fe XXIII 1.129 keV 
line, 33% for the Si XIII 1.865 keV line, 19% for the S XV 
2.460 keV line, 73% (!) for the Fe XXI 1.009 keV line; is only 



0 10 20 30 40 50 60 

Projected radius (kpc) 


Figure 13. Same as Fig. [n], but for the M87 cooling flow. 


15% for the Fe XXIV 1.17 keV and as small as 2% for the 
Fe XXV 6.7 keV line. 

Given the large number of resonant lines with a signif¬ 
icant optical depth from different elements and ionization 
species together with the proximity to us, the M87/Virgo 
cluster appears to be a perfect target for future X-ray po- 
larimetric observations. It is clear that the projected radial 
profiles shown in Fig. |l3| should be very sensitive to the tem¬ 
perature, density and abundance physical radial profiles. 


4 CONCLUSIONS 

• We have shown that the expected degree of polarization 
for the brightest X-ray lines from clusters of galaxies is of 
the order of 10%. 

• The degree of polarization is a function of projected 
distance from the cluster center, and the polarization vector 
is perpendicular to the projected radius vector. This implies 
in particular that the integrated X-ray flux from the whole 
cluster will be unpolarized (provided the cluster if spheri¬ 
cally symmetric). 

• It has also been demonstrated that for regular clusters, 
whose gas density distribution obeys the beta-law (jlj), the 
largest signal-to-noise ratio in measuring the polarization 
signal is attainable in the range of projected radii between 
0.5 and 10 core-radii. The degree of polarization is zero in 
the direction of the cluster center, while the photon flux 
from the outer parts is small making it difficult to measure 
the polarization. 

• Rich regular clusters as well as clusters with a dominant 
cooling flow are the best targets for future polarimetric ob¬ 
servations. 
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Why is it important to observe the polarization effect? 
First of all, this effect is unique in the sense that there is no 
other obvious mechanism that could lead to polarization in 
X-ray lines. Thus, a detection of polarized diffuse emission in 
an X-ray line would be a solid proof that resonant scattering 
is indeed in operation. 

As mentioned before, resonant scattering affects the 
“apparent” (inferred under the assumption of optically thin 
emission) element abundance radial profiles of clusters (Gil- 


fanov e; al. 1987a 


Shigeyama 1998). As a consequence, in 


order to determine the true abundance profiles from spec¬ 
troscopic measurements, one should be able to accurately 
subtract the distortions caused by resonant scattering, which 
can well be larger in magnitude than the true underlying ra¬ 
dial trend. This is a difficult task, because the surface bright¬ 
ness is an integral characteristic of a given line of sight which 
passes through regions with different densities and tempera¬ 
tures. One way to proceed is to compare the surface bright¬ 
ness profiles of at least two spectral lines corresponding to a 
given element (e.g. Ka and K[5 lines of iron, see Akimoto et 
al. 1997; Molendi 1998). Polarimetric measurements could 
significantly simplify such an analysis, for measuring of the 
polarized flux in a single line already allows one to estimate 
the contribution of resonant scattering to apparent abun¬ 
dance profiles. Apart from element abundances, the resonant 
scattering effects in principle can provide information on the 
velocity held within the cluster, because the optical depth in 
a given line is a function of the characteristic velocity of the 
corresponding ions. Lines of iron are particularly sensitive 
to the velocity of turbulent motions. 

Finally, the polarization effect can also be used for cos¬ 
mological purposes. Indeed, since the surface brightness of 
resonantly scattered emission and of its polarized compo¬ 
nent is roughly proportional to the optical depth along a 
given line of sight while the surface brightness of contin¬ 
uum thermal emission is proportional to the integral of the 
density squared along the line of sight, it is possible to deter¬ 
mine the distance to the cluster and consequently the Hub¬ 
ble constant ana logously as is done using resonance X- ray 
absorption lines ( Krolik fe Raymond 1988 ; Sarazin 198£ ) or 
the cosmic microwave background (the Sunayev-Zel’dovich 
effect). 

There are obvious requirements to a future X-ray po- 
larimeter. First, it should have the sufficient angular resolu¬ 
tion to probe individual parts of a cluster (we recall that the 
polarization vanishes when integrated over the whole clus¬ 
ter). Second, the CCD-type energy resolution is desirable 
in order to avoid significant contamination of the polarized 
emission of a resonance line by the unpolarized emission of 
neighbouring lines and the continuum. Finally, a large effec¬ 
tive area of the telescope is needed, larger than is required 
for the detection of unpolarized resonantly scattered emis¬ 
sion, since the expected degree of polarization is ~ 10%. 
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